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Abstract 



We analyze the asymptotic properties of an Euclidean optimization problem on the plane. 
Specifically, we consider a network with 3 bins and n objects spatially uniformly distributed, 
each object being allocated to a bin at a cost depending on its position. Two allocations are 
considered: the allocation minimizing the bin loads and the allocation allocating each object to 
its less costly bin. We analyze the asymptotic properties of these allocations as the number of 
objects grows to infinity. Using the symmetries of the problem, we derive a law of large numbers, 
a central limit theorem and a large deviation principle for both loads with explicit expressions. 
In particular, we prove that the two allocations satisfy the same law of large numbers, but they 
do not have the same asymptotic fluctuations and rate functions. 

Keywords: Calculus of variations; Central limit theorem; Euclidean optimization; Large deviations; 
Law of the large numbers; Wireless networks. 

1 Introduction 

In this paper we take an interest in an Euclidean optimization problem on the plane. For ease 
of notation, we shall identify the plane with the set of complex numbers C. Set A = 2(3\/3)~^/^, 
i = 1/— T (the complex unit), j = e^*'^/^ and consider the triangle T C C with vertices B2 = Xi, 
Bi = fB2, and S3 = jB2. Note that T is an equilateral triangle with side length A\/3 and unit 
area. We label by {1, . . . , n} n objects located in the interior of T and denote by Xj., k = 1, . . . ,n, the 
location of the A;-th object. We assume that {-^fc}fc=i,...,n are independent random variables (r.v.'s) 
with uniform distribution on T. Suppose that there are three bins located at each of the vertices 
of T and that each object has to be allocated to a bin. The cost of an allocation is described 
by a measurable function c : T ^ [0, cxd) such that ||c||oo := sup^g-j- c(x) < 00. More precisely, 
c{x) = ci{x) denotes the cost to allocate an object at a; G T to the bin in Bi; the cost to allocate 
an object at x S T to the bin in B2 is C2{x) = c{j'^x); the cost to allocate an object at x G T to the 
bin in B^ is c^^x) = c{jx). Let 



be the set of allocation matrices: if a^; = 1 the k-th object is affiliated to the bin in Bi. We consider 
the load relative to the allocation matrix A = {aki)i<k<n,i<i<3 £ -^n- 
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-^n — {A — {akl)l<k<n,l<l<3 '■ O-kl e {0, 1}, ttkl + ak2 + flfcs — 1} 




and the minimal load 



Pn = min Pn{A). 



B2 




Figure 1: The triangle T, the three bins and the n objects. 

Throughout this paper we refer to p„ as the optimal load. This simple instance of Euclidean 
optimization problem has potential applications in operations research and wireless communication 
networks. Consider three processors running in parallel and sharing a pool of tasks {l,...,n} 
located, respectively, at {^i, • • • ,^n} ^ T. Suppose that q(x) is the time requested by the l-ih. 
processor to process a job located at x G T. Then p„ is the minimal time requested to process all 
jobs. For example, a natural choice for the cost function is c{x) = 2\x — i.e. the time of a 
round-trip from Bi to x at unit speed. In a wireless communication scenario, the bins are base 
stations and the objects are users located at {^1, • • • ,Xn} ^ T. For the base station located at 
Bi, the time needed to send one bit of information to a user located at x G T is q(x). In this 
context Pn is the minimal time requested to send one bit of information to each user and 1/ Pn is 
the maximal throughput that can be achieved. We have chosen a triangle T because it is contained 
in the hexagonal grid, which is a good model for cellular wireless networks. 
For 1 < / < 3, we define the Voronoi cell associated to the bin at Bi by 

Ti = {xeT:\x-Bi\= min |x - 

l<m<3 

where Di = {ijt : t < 0} and, for / = 2, 3, A = {ift : t < 0}. Note that Ti U T2 U T3 = T and 
Ti n T2 = Ti n T3 = T2 n T3 = 0, i. e. {Ti, T2, T3} is a partition of T. Note also that G Ti. 

Throughout the paper, we denote by | • | the Euclidean norm on C, by i the Lebesgue measure 
on C and hy x ■ z the usual scalar product on C, i. e. x ■ z = 5R(x)3f?(z) + Q{x)Q{z). We suppose 
that the value of the cost function is related to the distance of a point from a bin as follows: 

For all X G T and / = 2,3, if |x — Bi\ < |x — Bi\ then ci(x) < q(x). (1) 

For example, if c(x) = /(|x — Bi\) and / : [0, oo) [0, oo) is increasing then ([1]) is satisfied. 
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In this paper, as n goes to infinity, we study the properties of an ahocation which reahzes the 
optimal load pn, and, as a benchmark, we compare it with the suboptimal load = where 
^ = i'(iki)i<k<n,i<i<3 is the random matrix obtained affiliating each object to its less costly bin: 

au = l{Xk G Ti). 

We shall prove that, using the strong symmetries of the system, it is possible to perform a fine 
analysis of the asymptotic optimal load. It will turn out that a law of large number can be deduced 
for the optimal and suboptimal load. More precisely, setting 

7 = / c(x) dx, 

we have 

Theorem 1.1 Assume ([7]). Then, almost surely (a.s.), 

lim — = lim — = 7. 

As a consequence, at the first order, the optimal and the suboptimal load perform similarly. 

The next result shows that, at the second order, the two loads differ significantly. We first 
introduce an extra symmetry assumption on c, namely, its symmetry with respect to the straight 
line determined by the points and Bi. If x = te*^ G T, t > 0, G [0, 27r], then its reffection with 
respect to the straight line determined by the points and Bi is te~*^~*3 G T. Formally, we assume 

c(te*^) = c(te"*^~*t) for ah 9 G [0, 2-k] and t > such that te*^ G T, 

and c is Lipschitz in a neighborhood of Di U D3. (2) 

Setting 

= c^{x) dx, 

we have: 

Theorem 1.2 Assume ([7]) and (0). Then, in distribution, as n goes to infinity, 

n-^l\pn-in)^G 

where G is a Gaussian r. v. with zero mean and variance — 7^. Moroever, in distribution, as 
n goes to infinity, 

i^^'^^^iPn - in) max{Gi, G2, G3} 

and 

n'^^^ipn - Pn) max{Gi, G2, G3} - -(Gi + G2 + Gs), 
where Gi, G2 and G3 are independent Gaussian r. v.'s with zero mean and variance o"^. Finally 

E[pn] = 727 + o{\/n) and E[p„] = 727 + m^/n + o{\/n), 
where m = E[max{Gi, G2, G3}] > depends linearly on a. 
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Theorem 11.11 states that p„ is asymptotically optimal at scale n, but Theorem 11.21 says that it is 
not asymptotically optimal at scale \/n. In the proof of Theorem 11.21 we shall exhibit a suboptimal 
allocation which is asymptotically optimal at scale \/n (see Proposition 13. ip . 

We shall also prove a large deviation principle (LDP) for both the sequences {pn/"'}n>i and 
{p„/n}„>i. Recall that a family of probability measures {/in}n>i on a topological space (M, Ta/) 
satisfies a LDP with rate function / if / : M — > [0, oo] is a lower semi-continuous function such that 
the following inequalities hold for every Borel set B: 

— inf /(y) < liminf — log^„(-B) < limsup — log /U„(-B) < — inf_/(y), 



n^oo n n^oo n y^B 



where B denotes the interior of B and B denotes the closure of B. Similarly, we say that a 
family of M-valued random variables {ln}n>i satisfies an LDP if {nn}n>i satisfies an LDP and 
/i„,(-) = PiVn G •)■ We point out that the lower semi-continuity of I means that its level sets 
{y G M : I[y) < a} are closed for all a > 0; when the level sets are compact the rate function /(•) 
is said to be good. For more insight into the large deviations theory, see, for instance, the book by 
Dembo and Zeitouni [4j. 

We introduce an assumption on the level sets of the cost function: 

^(c~^({t})) = forallt>0, (3) 

an assumption on the regularity of c: 

c is continuous on T, (4) 

and two further geometric conditions: 

c(Si) < c(x) < c(0), for any x G Ti \ {0, Bi}, (5) 



Cl(x)c2(x)c3(x) c(0) /" . s , , r 1 

, , , ^ /^ /X , ^ < ^ < / c z dz, for any X G T \ {0}. 

Ci{x)C2{x) + Ci{x)C3{x) + C2{X)C3{X) 3 Jj 



Assumption ([5]) fixes the extrema of the cost function on Ti. The left hand side inequality of 
([6]) imposes that is the most costly position in terms of load (for a more precise statement, we 
postpone to ([37]) ). For G M, define the functions 

A{9)=log(^3j e^^^'^dx^ and A{9) = log(^j e^"(^') dx + 2/3^ 

and, for y G M, their Fenchel-Legendre transforms 

A*(y) = sup(%-A(6')) and A*(y) = sup(% - A(0)). 
em 6»eK 

The following LDPs hold: 

Theorem 1.3 Assume 0, ® and ®. Then 

(*) {Pn/n}n>i satisfies an LDP on R with good rate function 

J(y) = I '^y^ (c(i?i)/3,c(0)/3) 

I +00 otherwise. 

(ii) {7?„/^^}n>l satisfies an LDP on M with good rate function 

A*(3y) if y G (c(i?i)/3,7] 
Jiy)={ A*(y) ifyG(7,c(0)) (8) 
-|-oo otherwise. 
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Next proposition gives a more explicit expression for the rate functions. 

Proposition 1.4 Assume @ and c continuous at and Bi. Then A* and A* are continuous 
on (c(-Bi), c(0)) and 

Av,l = / y^y-^(^y) if c(i?i) < y < c(0) 

\ +00 if c{Bi) >yoiy> c(0) 

where 9y is the unique solution of 



/^^ c(x)e'''=(^) dx 
/j,^ e^'=(^) dx 

(n) 



2/- (9) 



A (y) 



yr]y - A{i]y) if c{Bi) <y< c(0) 
+00 if c{Bi) > y or y > c(0) 

where r]y is the unique solution of 

c(x)e^'=(^) dx 



/j,^ e^'=(^) dx + 2/3 



(10) 



(m) Ifj<y< c(0)/3 i/ien A*(y) < A*(3y). 

Note that: J{y) = A*(3y) except possibly at y G {c(i?i), c(0)}; J(?/) = A*(3y) on {—00,^] except 
possibly at y = c(i?i), and J{y) = A (y) on (7,00) except possibly at y = c(0). These gaps are 
treated in Proposition 14.41 with extra regularity assumptions on c. See Figure [2] for a schematic plot 
of the rate functions. A simple consequence of Theorem 11.31 and Proposition 11.41 is the following: 

n^<^ log P{p^>nt) J{t) n^oo P{p^>nt) ' ; 

In words, it means that the probability of an exceptionally large optimal load is significantly lower 
than the probability of an exceptionally large suboptimal load; although, on a logarithmic scale, the 
probability of an exceptionally small optimal load does not differ significantly on the probability 
of an exceptionally small suboptimal load. It is not in the scope of this paper to discuss the trade- 
off between algorithmic complexity and asymptotic performance. Moreover, we do not know if 
the allocation that is asymptotically optimal at scale ^/n used in the proof of Theorem 11.21 (see 
Proposition 13. ip has the same rate function than pn/n. 

Unlike it may appear, we shall not prove Theorem 1 1.3 1 by first computing the Laplace transform 
of pn and p„ and then applying Gartner-Ellis theorem (see e. g. Theorem 2.3.6 in [3]). We shall 
follow another route. First, we combine Sanov theorem (see e. g. Theorem 6.2.10 in [3]) and the 
Contraction Principle (see e. g. Theorem 4.2.1 in [4]) to prove that the sequences {pn/n}n>i and 
{p„/n}n>i obey a LDP, with rate functions given in variational form. Then, we provide the explicit 
expression of the rate functions solving the related variational problems. It is worthwhile to remark 
that, using Theorem 1 1 . 3 1 and Varadhan lemma (see e. g. Theorem 4.3.1 in [4J) it is easily seen that 

lim ilogE[e^^"] = J*(0) and lim i log £[6^^^"] = J*(6'), V6l G M 

n—fOO n n— >oo fi 
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where J* and J are the Fenchel-Legendre transforms of J and J, respectively. A nice consequence 
of Theorems 1 1 . 1 1 and 1 1 . 2 1 is that, in terms of law of the large numbers and central limit theorem, p„ 
has the same asymptotic behavior as 



^ 3 n 

-Y,Y.^{XkeTi]ci{Xk). 



Pn 2 

1=1 k=l 

Moreover, if the cost function satisfies extra regularity assumptions (see Proposition I4.4p . by Theo- 
rem [L3] and the Gartner-Ellis theorem, we have that p„ and pn have the same asymptotic behavior 
even in terms of large deviations. 

As it can be seen from the proofs, if the left hand side of assumption ^ does not hold then we 
have an explicit rate function J{y) only for y < c(0) /3. If neither the right hand side of assumption 
([6]) holds, then we have an explicit rate function J{y) only for y < yo for some yo > 7. We also 
point out that the statements of Theorems 11.2111.31 concerning do not require the use of ([2]) and 
©. 

In wireless communication, the typical cost function is the inverse of signal to noise plus inter- 
ference ratio (see e.g. Chapter IV in Tse and Viswanath [9]), which has the following shape: 

, , a + min{6, |x - i?2|~"} + niin{6, |x - i?3|"°} ^ 
^ ^ mm{6, \x - 

where a > 2, a > and 6 > (A^/3/2)-° (recall that A = 2(3\/3)-^/2 ^nd A\/3 = \Bi - ^al). We 
shall check in the Appendix that this cost function satisfies ([T|), ([3]), ^ and ([5]). Moreover, 
the first inequality in ([6]) will be checked numerically and, for arbitrarily fixed a > 2 and a > 0, 
we shall determine values of the parameter b > (A\/3/2)~" such that the second inequality in ([6]) 
holds. 

The remainder of the paper is organized as follows. In Section [2] we analyze the sample path 
properties of the optimal allocation and we prove Theorem 11.11 In Section [3] we show Theorem 
II. 2i Section U] is devoted to the proof of Theorem 11.31 and Proposition 11.41 In Section \5\ we discuss 
some generalizations of the model. We include also an Appendix where we prove some technical 
lemmas and provide an illustrative example. 



2 Sample Path Properties 

2.1 Structural properties of the optimal allocation 

Throughout this paper we denote by Mb(T) the space of Borel measures on T with total mass less 
than or equal to 1 and by Mi(T) the space of probability measures on T. These spaces are both 
equipped with the topology of weak convergence (see e. g. Billingsley [1]). For a Borel function h 
and a Borel measure jj, on T, we set n{h) = Jj h{x) iJ,{dx). Consider the functional from Mfe(T)'^ to 
M defined by 

(j){ai,a2,a3) = max(ai(ci), a2(c2), 03(03)). (11) 

Letting a\Q denote the restriction of a measure a to a Borel set B, we define the functionals $ and 
^ from Mi(T) to M by 

$(a) = inf (/>(«!, 0:2, as), 

and 

*(a) = (/)(aiTri,a|T2>a|T3)- 



6 



J 



J 




3 



7 £{0) 
' 3 



c(0) 



Figure 2: The rate functions J and J. 



Note that if 6x denotes the Dirac measure with total mass at x S T, then 

— / 1 " ^ 



n 



n 



(12) 



k=l 



Lemma 2.1 Under assumption Q we /laue i/iai (/) is continuous on M{,(T)^ and ^' and $ are 
continuous on Mi(T) (/or i/ie topology of the weak convergence) . 

The proof of Lemma 12.11 is postponed in Appendix; the continuity of (j) and ^ is essentially trivial, 
the continuity of $ requires more work. Define the set of matrices 



and 



'^n = {B = {bkl)l<k<n,l<l<3 ■ bkl G [0, 1], 6fcl + bk2 + bk3 = 1} 



minp„(S). 



Given a matrix B = [bki] G Tin, we define the associated measures (ai,a2,Q!3) G M;,(T)^ by setting 
ai = (1/n) 'Yllt=i bki^Xu i} = li 2, 3). Due to this correspondence, it is straightforward to check that 



(13) 



k=l 



Next lemma is a collection of elementary statements, whose proofs are given in Appendix. 

Lemma 2.2 Fix n > 1 and let B* = (b^li) G "Bn be an optimal allocation matrix for pn- Then: 

(i) For all a G Mi(T), there exists (ai, 02,03) G Mfc(T)^ such that a = 01 + 02 + 03 and <I>(o) = 
4>{ai, 02, 03). Moreover, whenever such equality holds, we have oi(ci) = 02(02) = 03(03) and, 
in particular, 



^biMXk) = Ebi,c2iXk) = Y,KMXk). 



k=l 



k=l 



k=l 
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(a) If assumption holds then 

Pn - 3\\c\\oo < Pn < Pn a.S.. 

{Hi) If assumption (0) holds then the sequences {pn/n} and {pn/n} are exponentially equivalent. 
For the definition of exponential equivalence see [4j. 

2.2 Proof of Theorem [TT] 

The law of large numbers yields, for all / = 1, 2, 3, 



n 

lim -Yci{Xk)l{Xk£Ti} = / ci{x)dx 



7 a.s. 



Therefore from the identity 



^ = max - G TJ, 

n KK3 n ^ 



n 



k=l 



we get lim„_^oo Pn/^ = 7 ^--S- We also have to prove that lim„^oo Pn/n = 7 a-s.. Let A = (aki) G An 
be an allocation matrix. By assumption ([T]), if x G T/ then q(x) = mini<m<3 Cm{x). Therefore 



3 n 

Spn{A) >Y.Y1 ^kici{Xk) >Y.Y1 ^'(^'^^ ^ 3 i^^3 ( E G j . (14) 

/=1 k=l 1=1 XfeGT, \fc=l / 



So taking the minimum over all the allocation matrices we deduce: 

1^^3 (j2ci{Xk)l{Xk G Ti}^ < Pn < Pn- 

Thus by applying the law of large numbers, we have a.s. 

7 < lim inf — < lim sup — < 7. 

n^oo n n^oo n 

Remark 2.3 Assume that conditions ([T|), ([3]) and ^ hold. By Theorem 1 1.1 1 we have lim^^oo T^n/^ — 
7 a.s.. So by Lemma |2.H equation (fT2]) and the a.s. weak convergence of (1/^^) l^fc=i (^x^ to 
£ we get ^'(^) = 7. Similarly, using equation (fT3l) in place of equation (fT2l) . we deduce that 
lim^^oo /On/"- = ^{£) a.s.. So by Lemma \2.2{ ii) limn^oo Pn/n- = ^J'l^) a.s., and by Theorem 11.11 
we have ^{£) = 7. 

3 Proof of Theorem D 

Consider the random signed measure 

1 " 

Wn = Vn {pn - £) where /i„ = - • 

k=l 
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The standard Brownian bridge Vl^ on T is a random signed measure specified by the centered 
Gaussian process {W{f)} (indexed on the set of square integrable functions on T, with respect to 
i), with covariance given by 

E[W{f)W{g)]=i{fg)-£{mg), 
see e.g. Dudley [5]. By construction 

'Pn = n max ( / ci{x)iin{Ax, 
KK3 \ T, 



or equivalently 



/Jl l<i<3 \Jj 

Let / be a square integrable function on T, in distribution, as n — > oo 

\ n 



max / ci{x)Wn{dx) . (15) 



indeed by the central limit theorem Wn{f) converges in distribution to a Gaussian r.v. with zero 
mean and variance equal to i{f'^) — ^^(/), which is exactly the law of W{f). Using Levy continuity 
theorem and the inversion theorem, we have, in distribution, for all square integrable functions /i, 
/2 and /s: 

{Wn{fl),Wn{f2),Wn{h)) ^ {W {h) ,W {f2) ,W [h)) . 

Therefore, by the continuous mapping theorem we have, in distribution, as n goes to infinity, 

'Pn-n^ - max f / ci{x)W{dx) \ . (16) 



n KK3 



We shall show later on that the r.v. in the right-hand side of (jl6p has the claimed distribution. Now 
we consider the optimal load /?„. By the second inequality in (jl4p we have 



3 

1=1 J'^l 



and therefore 

^Pn^ [ ci{x)Wn{dx). (17) 

The following proposition is the heart of the proof. It will be showed later on. 

Proposition 3.1 Under the assumptions of Theorem M.iA for any 1/4 < a < 1/2, there exists 
an allocation matrix A = {aki)i<k<n,i<l<3 £ -An with associated load pn = Pni^) such that with 
probability at least 1 — Li exp(— Lgn-*^"^"), 



for some positive constants Lq and Li, not dependent of n. 
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Using this result, /5„ > /?„ and p7|) . we have that with probability at least 1 — Li exp(— Lgn 



Pn - n-f 



n 



[ Ci{x)Wn{dx] 

l</<3 



(18) 



Therefore, as n goes to infinity, in distribution 

l</<3 



n 



I yZ [ Ci{x)Wn{dx) ^ 0. 

"J ^"^^ It, 



The continuous mapping theorem yields 



So combining these latter two limits we get, as n goes to infinity, 

3 



Pn-nj 1 



n 



Y [ ci{x)Widx), 
1=1 '''^i 



i. e. n~^^'^{pn — nj) converges weakly to a centered Gaussian random variable with variance a'^/3 — 
7^. We have considered so far, the normalized sequences pn and p„ separately. However, we can 
carry the same analysis on the normalized difference p„ — /)„. More precisely, by equation (fT5l) we 
have a.s. 



Pn - Pn 



< 



n 



Pn - 



max 



max 



n 



ci{x)Wn{dx) 
ci{x)Wn{dx 



1<Z<3 



Y / Cl{x)Wn{dx) 



+ 



Pn-nj 1 



n 



KK3 



Y / Ci{x)Wn{dx) 



Pn - _ 1 
3 



n 



Y [ Cl{x)Wn{dx) 
KK3 



Thus, by equation ([TH]) . we obtain, with probability at least 1 — Li exp(— Lon 



Pn - Pn 



n 



max I / ci{x)Wn{dx 
1<K3 VA, 



^ Y I Ci{x)Wn{dx) 
KK3' 



- 3 



Therefore, in distribution, as n ^ oo, 



Pn - Pn 



n 



max I / ci{x)Wn{dx] 
l<i<3 V./t, 



\Y I ^li^Wn{d2 



0. 



The continuous mapping theorem yields 

Ci{x)Wnidx) 



max 

l</<3 



y cKx)iy„(dx))-l Y i 



max 



q(x)Ty(dx)V^ Y [ Ci{x)W{dx) 
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and therefore, in distribution, as n oo, 

max ( / ci{x)Widx)) -\Y. I ci{x)W{dx). 

For / G {1,2,3}, set 

3 

iV/= / ci[x)W(Ax)-\y^l ci[x)W(Ax). 
Hi -3 ^ Hi 

By definition is a centered Gaussian process indexed on the set of square integrable func- 

tions, therefore N = (Ni, N2, N^) follows a multivariate Gaussian distribution with mean 0. A 
simple computation shows that the covariance matrix of A'^ is 




It implies that N has the same distribution as 

(Gi - (Gi + G2 + G3)/3, G2 - (Gi + G2 + G3)/3, G3 - (Gi + G2 + G3)/3) 

where Gi,G2 and G3 are independent Gaussian r. v.' s with mean and variance o"^. 

It remains to compute the asymptotic behavior of the expectation of the loads. A direct com- 
putation gives, for any / = 1, 2, 3, 



E 



Ci{x)Wnidx] 



3 9n ~ 3 ' 



Thus the sequences {jj^ q(x)VF„(dx)} (Z = 1,2,3) are uniformly integrable. This implies that the 
sequence |maxi<K3 ^ jj^ ci{x)Wn{dx)^ | is uniformly integrable and so using equation (fTSi) we have 



lim E — nj] / \fn = lim E 

n— >oo n- — 

= E 



max I / ci{x)Wn{dx\ 

1<K3 \Hi 



max 

KK3 



ci{x)W(dx) 



m. 



Now we give the asymptotic behavior of E[/3„]. Note that by (jlSp we have 



E 



V / Q(x)iy„(dx) 

■TTt-n H, 



1<Z<3 



, Pn - nj 



n 



V / Cl{x)Wn{dx) 



1{| . . . I >n 



l/2-2a 



< n^/^-2« + lO||c||ooi>iVnexp(-Lon^~2a) 

= n^/2-2a ^ Li V^exp(-Lon^-2") (19) 

where the latter inequality follows since 7 < ||c||oo, Pn < ||c||oof^ and | Jrj,^ q(x)IV„(dx)| < 2||c||oo-v/^- 
Clearly, the term in ([T9|) goes to zero as n — > 00. Therefore, since E ci{x)Wn{dx) = 0, we have 

lim E[p„ — n'^]/y/n = 0. 
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Proof of Proposition COl We start describing the allocation matrix A. For l,m G {1,2,3} and 
t G [— A\/3/2, A\/3/2], denote by Bim{t) the point on the segment BiB^ at distance t + \^/^/2 from 
El. We extend the definition of Bim{t) for all t G [— A\/3, A\/3] by following the edges of T. More 
precisely, we set 

53i(A\/3 + t) if t e [-A^/3, -A\/3/2] 



R (f\ = \ B-ii{\V?, + t) if t e [-AV3, -A^/3 
^'^ ^ lB23(A^/3-t) if t G [AV3/2,A^/3]. 



For l,m £ {1,2,3}, Bifn{t) is defined similarly by a circular permutation of the indices. For t 
{t\t^,t^) e [-XV3,XV3f, let 

C\t) = {0} U ({z G C : z • (Si2(i^)e-*"/2) > 0} n {z G C : z • {B3i{t^)e'^/^) > 0} 

be the (possibly empty) cone delimited by the straight straight line determined by the points and 
B3i{t^). We define ^^(t) = C^{t)nT. Similarly, let D'^{t) = C'^{t)nT and D^{t) = C3(t)nT with 

C^t) ={zeC:z- (5i2(t^)e*"/2) > 0} n {z e C : z ■ {B2^{t^)e-'^'^) > 0}, 
C^{t) = {z£C:z- (B2z{t^y^'^) > 0} n {z G C : z • {B^i{t^)Q-'^/^) > 0}. 

By construction, the sets D^{t), D'^(t) and D^{t) are disjoint and their union is T. For / G {1, 2, 3}, 
set 

n 

pU^) = ^ci{x,)i{XkeD\t)} 

k=l 

and consider the following recursion. At step 0: for to = (0,0,0), define 

mo = argminp^(to) 

l</<3 

(breaking ties with the lexicographic order) and 

Mo = argmaxp^(to) 

l</<3 

(again breaking ties with the lexicographic order). If /3^^°(to) — /o™°(to) < 2||c||oo, the recur- 
sion stops. Otherwise, p^^°(to) — pJ^°(to) > 2||c||oo and there is at least one point Xi {i = 
l,...,n) in D*^''(to). Note also that, a.s., for all 9 G [0,27r], there is at most one point of 
{Xi,--- ,Xn\ on the straight line (xe*^,x > 0). As a consequence there exists a random vari- 
able < ti < A-v/3 such that, a.s., there is exactly one point Xi {i = l,...,n) in the triangle 
with vertices {0, -BmoMo(*i)) -BmoMo(O)} for < ti < Aa/3/2, or in the polygon with vertices 
{0,S^oMo(ii),5Mo,5moMo(0)} for AV3/2 < ti < \^/^. We then set ti = (ti,0,0) if mo = 1, 
Mo = 2; ti = (-ti,0,0) if mo = 2, Mo = 1; ti = (0,ti,0) if mo = 2, Mo = 3; ti = (0, -ti,0) 
if mo = 3, Mo = 2; ti = (0,0, -ti) if mo = 1, Mo = 3; ti = (0,0, ti) if mo = 3, Mo = 1. By 
construction, we have 

Pr(ti) < Pf°(ti) , max^pL(ti) < max pi(to) and mm^p^(ti) > mm^p^(to). 
At step 1: define 

mi = argminp^(ti) 

1<«<3 

(breaking ties with the lexicographic order) and 

Ml = argmaxpj^(ti) 

KK3 
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(again breaking ties with the lexicographic order). Similarly to step 0, if p*^^ (ti) — (ti) > 2||c||oo, 
then there is at least one point of {Xi, • • • , X„} in D^^^ (ti) and we build the random vector t2. The 
recursion stops at the first step k > such that 

p^'Htk)-pT(.tk)<2\\c\\^, 

(where m^, and are defined similarly to niQ, mi, . . ., Mq, Mi, . . . and ti, t2, . . .). As we shall 
check soon, the recursion stops after at most n steps. When the recursion stops, say at step kn < n, 
we set = D^t^^) and t„ = t^^. The allocation matrix A is defined by allocating to the bin 
in Bi if Xfc G Dl, i. e. 

A = iaki)i<k<n,i<i<3 where aki = l{Xk G -Oi}. 
By construction, we have for all l,m £ {1,2,3} 

|pL(tn)-p;r(tn)| <2||C|U. (20) 

We now analyze the recursion more closely. Assume that at step we have mo = 3 and Mq = 1, 
i. e. pl{to) > /o^(to) > Pni^o)- Then, for ah k < kn, 

/Oi^(tfc) > /O^(tfc) - ||c||oo and /9n(tA:) < /O^(tfc) + ||c||oo- (21) 

Indeed, if for all k < kn, ruk = 3 and = 1, there is nothing to prove since |p^(tfc+i) — Pn{tk)\ < 
||c||oo- Assume that there exists k < kn such that m^ 7^ 3 or 7^ 1. We define 

kQ = mm{k > 1 : m^. 7^ 3 or 7^ 1}. 

For concreteness, assume for example that M^g 7^ 1. By construction, feg — 1 < kn so that 
Pn(tfeo-i) > /'n(tfco-i) + 2||c||oo- Since p^(tfco_i) > pl{tko-i) > /)i^(tfc„_i), we deduce that M^,,, = 2 
and niko ~ Recall that, for k < kn, /5^*(tfc) — ||c||oo < Pn''{^k+i) < Pn'^i^k)- Thus, for k = k^ — l, 
from /3^(tfcJ < pl{tko) = Pn(tfco-i) < /'n(tfco-i)' we obtain 

pli^ko) - \\c\\oo < pi{-tko) 

Similarly, for k < kn, PnH^k) + ||c||oo > PnH^k+i) > Pn^i^k)- Thus, from pU^ko-i) < pli^ko) = 
/9^(tfcg_i), we have 

piitko) < l|c||oo +/?^(tfco)- 

We have proved so far that the inequalities in ([2T]) hold for all k < k^. Since |/9^(tfc+i) — Pn{tk)\ < 
||c||oo and pl(tko-i) - pl(tko-i) > 2||c||oo we get 

/On(tfco)-Pn(tfco) >0- 

Thus niko = 3 and /9^(tfcJ < ^^(tfcj < pl{tko)- Define 

ki = min{A;„, min{/c > fco : m^, 7^ 3 or 7^ 2}}. 

For k = kf), - ■ ■ ,ki — 1, Pn{tk+i) < Pni'^k) and /?^(tfc+i) = p^(tfc) is constant, so the left hand side 
inequality of (f2TI) holds. Also, since /ci <kn,iork£ {kQ + 1,--- , ki — 1}, Pn{tk) < /9^(tfc)+4||c||cxD. So 
finally, (pT|) holds for /c = 0, • • • ,ki. Moreover, if ki < kn, then M^^ = 1 and m^^ = 3. Indeed, as 
above, p^(tfc,_i) - /)^(tA,.j__i) > 2||c||oo implies 

^n(tfci) > Pn(tfci)- 
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So Mfc^ 7^ 3 and m^^ 7^ 2. If = 1 and Mj^^ = 2, then we write, by (j21|) . 

/3n(tfci) + l|c||oo > pi{tki) > /3n(tfci) > /5n(tfcl)- 

So ki = kn, a contradiction. Therefore, we necessarily have M^^ = 1 and = 3. By recursion, 
it shows that for all k < kn, = 3. Hence, at each step one point is added to the bin at B^. No 
point is added to the bins at Bi and B2, points may only be removed from the bins at Bi and 
B2. Since there are at most n points, we deduce kn < n, as claimed. Also, since -D'(to) = T/, we 
obtain, for all A; = 1, • • • , kn, T3 C -D^(tfc), T2 5 D^{tk) and Ti D D^(tk). The other case, where 
rrifco = 2 could be treated similarly. So more generally, if, at some step, / = nik then / 7^ Mj for all 
k < j < kn, and conversely, if / = Mk then / 7^ for all k < j < kn- It implies that D\tk) is a 
monotone sequence in k. Since D'(to) = T/, for all I € {1, 2, 3}, 

Dl<^Ti or T^CZ)^. (22) 

Assume now, that t\ > zn~'^ with z > then, from ([22]), Ti C D\ and Z)^ c T2. For t G M, define 
the set V^{t) = 0, 0)\Ti. On the event > -zra""} we have 

Pi(t„) > n / c(.).„(d.) + „ / c(.„„(d., a.d A(K) < n f c.(.),„(d.). 

So, by inequality (pOj) . we deduce that on {t^ > zn~°'} 

I c(x)//„(dx) + / c(x)/i„(dx) < / C2(x)/i„(da;) + ^i^^^. 

Or equivalently, 

{t^ > zn-°} C J ^ / c(2;);u„(dx) < / C2(x)W„(dx) - / c(x)Ty„(dx) + ^M^i I . 

[ Jv^{zn-") JT2 Jfr V" J 

(23) 

Let ^ be a Borel set in T, by Hoeffding concentration inequality (see e. g. ^) we have, for all s > 
and / e {1,2,3}, 

P (^j ci{x)ixn{<^x) - j^ci{x)dx >s^< exp{-Kos^n), (24) 
P ci{x)ij,nidx) — J ci{x) dx < —s^ < exp(— ZCos^n) (25) 
where Kq = 2||c||j^^. Taking s = yn~°', where y > 0, we have 

P {^j^ ci{x)Wn{dx) > < expi-Koy^n^'^''), (26) 

P ci{x)Wn{dx) < -yn^-°^ < exp{-Koy^n'-^''). 
Similarly, by ([25]) we deduce, for s > 0, 

/ c{x)fin{dx) < / c{x)dx — s\ < exp{—Kos'^n). 



P 
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By assumption ([T]), there exists cq > such that c(x) > cq, for all x S V^{zn If < s < A-v/3/2, 
the area of V^{s) is equal to As/4. Therefore, for all < z < A\/3n"/2, 

Kizn-" < [ c{x) dx < K2zn~'^, 

with Ki = coA/4 and K2 = ||c||ooA/4. So, taking s = Kizn~°^ /2, we get, for all < z < Av^^^"; 

PlV^f c{x)iJr^{dx) < ^zn^-A < exp(-i^3z2ni-2°) (27) 

where = KqKi /4. Similarly, for t > 0, define 

U\t) = (l)'(t,0,0)\T,) IJ (D'(-t,0,0)\T,(,)) , 
where o" = (1 2 3) is the cyclic permutation. By (p^ we have, for all s > 0, 

PI / c{x)nn{dx)> / c(x) dx + s < exp(— Kos^JT-)- 

\Jm(zn-") Jm(zn-") J 

Thus, setting s = zn~'^, we get 

P {^Xn{U\zn-'')) > Kizn-'') < exp(-Ko^^n^-2°) (28) 
with K4 = 1 + 2i^2- Now, note that by ()23p . from the union bound, for y > 0, 

{ti> zn-"}<z\^ [ c{x)fin{dx)<yn-2-A\\{- [ ci(x)PF„(dx) + ^ > ^yn^^ 

U C2(x)M^„(dx) + ^ > ^yn^-| . 
Now take y = Kiz/2, by ([26]) and ([271) we deduce, if 4||c||oon""i/^i"^ < z< A\/3n" 

> zn-") < exp (-Kgz^ni-^") + 2exp (^-:^ni-2a (^^^ _ 4||c||^n"-i)^^ 
< 3exp f-i^sn^-^a ^^^^ _ 4||c||^„"-i)2 



with = mm{K3K^ ^, Xq/IG}. Therefore, by symmetry, for all n and z > such that 4||c||oo?i-" ^ ^ 
z < AV^n"/2 



P (^max > zn^"^ < 18 exp (^-i^gn^"^" (i^iz - 4||c||oon""^)^) . 
Note that pn = Pn{A) = maxi<K3 /)J^(t^), so by (pOj) we have 

3pn - 4||c||oo < P^(tn) + P^(tn) + />n(tn) < ^Pn- 

Subtracting 3^/n'y, it follows 

3—^= ^ < Q(x)^„(dx) - 7j < 3^^. 
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(29) 



Then we subtract the quantity 

3 



J2 / ci{x)Wn{dx) = VnJ2 ( / Q(2;)^„(d; 
1=1 '''^i 1=1 ^'''^i 



x) -7 



and we get 



,Pn - n-f 



V / ci{x)Wn{d. 



1=1 



+ 



4||c|| 



n 



(30) 

Set Cmin(ic) = miii(ci (x) , C2 (x) , C3 (x)) , and note that if x G T/ then Cmin{x) = c;(x). If tj^ > 0, we 
set V}, = = \ Tz, and, if tj, < 0, we set = Dl^^^ \ Tu where a = (12 3) is the cydic 

permutation. So 

3 3/1 3 

/ Q(x)^„(dx) - / Cl{x)fln{dx) {Cl{x) - Cminix)) Hn{dx) 

1=1 ■'Dl Jf, Jdi 

= ^ > 0} j^^ (q(x) - Cmin(2;)) /U„(dx) + ^ < 0} j^^ (c<^(0(3;) " Cmin(x)) ^„(dx). 



l=\ 



(31) 



Note that if x G T^, with m ^ I, then |q(x) — Cmm(2;)| = |q(2;) — Cm{x)\. For example, assume 
/ = 1, m = 2 and X = te'f e T2, with < 6* < 7r/3, we then have 

|ci(x) - c^i„(x)| = |ci(x) - C2(x)| = |c(te*t+*^) - c(te*t+*^e-*^)| 

= |c(ie*t+*^) -c(te-^i+^^)|. 

By the symmetry assumption ([2]), we deduce 

|ci(x) - c^i„(x)| = |c(te^t+^^) - c(te^|-^^)|. 

Again by assumption ([2]), c is Lipschitz in a neighborhood of Di U D^. Letting L > denote the 
Lipschitz constant, if x is close enough to Di, say the distance d(x,Di) from x to Di is less than 
or equal to < e < A-v/3/2, we have 

|ci(x) - Crnm(2;) < Lt|e*t+*^ - e^i-^^l = Li|e^^ - e"^^! 

= 2Ltsm9 = 2Ld{x,Di). 

By symmetry, for all I £ {1, 2, 3}, if d(x, D^) < e, then 

\ciix) - Cmin(2;)| < 2Ld(x, A) and \c„(l){x) - Cmin(x)| < 2Ld(x, A)- 

Fix a S (1/4, 1/2), z > and choose n large enoug h so that 4||c||oon°"iii:f ^ < z < en". Then, by 
([29]) with probability at least 1— 18exp (^—K'^n^~'^" {Kiz — 4||c||oof^"^"'^)^^ , we have maxi<z<3 < 
zn~'^. On this event, if x G V''(iJi) then d(x, Di) < zn~" < e. It follows by (j3ip that, with probability 
at least 1 - 18exp (-K^n^-^"^ (Kiz - 4||c||oon"-i)^ 



/ Q(x)/i„(dx) - / Q(x)/u„(da 
z=i -^Di, 1=1 J'^i 



1=1 

3 

< 2Lzn5'°^/i„(;7'(zn-")), 
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By i^B, with probability at least l-3eyip{-Koz^n^-^'^), itholdsJ2Li f^niU^zn''^)) < 3Kizn-°'. Us- 
ing that for all events A,B it holds P{A D B) > 1 — P(A^) — P{B^), we obtain, for all n large 
enough so that 4||c||oo'T'"~^-f^r^ ^ z < en", 

with probability at least 1—21 exp ^—Kqu^^'^" {^Kiz — 4||c||oo''^"~^)^^ , where Kq = min{i^o-f^i^^, K^}. 
this latter inequality and ([30|) . with the same probability, 

< 12Li^2-z^ra^-2" + 4||c||oon"5. 

Fix z = {24:LK2)^^ so that 12Li^2-z^ = 1/2. Then there exists no such that, for all n > no, 
4||c||cxDn"~^i^['"'^ < z < en" and 8||c||oo?i- ^ < na"^". Then, for all n > no, 

< ni/2-2a ^32) 

with probability at least 

1 - 21 exp ^-ivTen^-^" (|Ki(24LA"2)"5 - 4||c||oo<~^) = 1 - -fsTyexp (-i^g"^"^") ■ 

Finally, we set Lq = Kg and Li = max{Kj, Kg}, where Kg = exp (iiTsnJ-^") . With this choice of 
Lq and Li, (I32p holds for all n > 1 with probability at least 1 — Li exp(— Lon^~^"). 



V / q(x)/i„(dx) - V / 

1=1 -^Di, Jt, 



n 



3 

V / Ci{x)Wn{dx) 



n 



3 

V / Q(x)VF„(dx) 



4 Large deviation principles 

In this section we provide LDPs for the optimal and suboptimal load. Letting <C denote absolute 
continuity between measures, we define by 



1^ +00 otherwise 



the relative entropy of G Mi(T) with respect to the Lebesgue measure i. Moreover, if / is a non- 
negative measurable function on T, we denote hy if the measure on T with density /. In particular, 
if fj f{x) dx = 1, we set 

Hif) = Hiif\e)= [ f{x)logfix)dx. 

4.1 Combining Sanov theorem and the Contraction Principle 

Next Theorem 14. II follows combining Sanov theorem and the Contraction Principle. 

Theorem 4.1 Assume (0) and 0. Then 

{i) {/5n/n}„>i satisfies an LDP on M with good rate function 

J(y)= inf H(a\i). (33) 
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{ii) {/0„/^^}n>l satisfies an LDP on M with good rate function 

J(y)= inf H(a\i). (34) 

ae3Vti(T):*(a)=j/ 

Proof. By Sanov theorem the sequence Y^'^=i <5xj}n>i satisfies an LDP on Mi(T), with good rate 
function H{- \ Recall that the space Mi(T), equipped with the topology of weak convergence, 
is a Hausdorff topological space (refer to [1]). By Lemma 12.11 the function $ is continuous on 
Mi(T). Therefore, using and the Contraction Principle we deduce that the sequence {pn/n'}n>i 
satisfies an LDP on M with good rate function given by (j33p . Consequently, by Lemma [2.2r iu) and 
Theorem 4.2.13 in [4], {pn/n'}n>i obeys the same LDP. The proof of (ii) is identical and follows 
from (ini). □ 



Remark 4.2 It is worthwhile noticing that one can prove Theorem 14.11 also applying Lemma |2.H 
Lemma l2.2( ni) and the results in O'Connell [7|. 

4.2 Computing A* and A* 

In this subsection we compute the Fenchel-Legendre transforms A* and A*. 
4.2.1 Proof of Proposition 11.41 

We only compute A* in (i). The expression of A* in (ii) can be computed similarly. Clearly, for 

^ " }^^e(><-)dx ' 

and 

f „ pdcix) ( f Jc{x) \ 
A"{d)= c\x)- —^dx-{ c{x)-. ^TT-^dx >0, 

(the strict inequality comes from the assumption that c(-) is not constant on Ti). Therefore, the 
function A' is strictly increasing. Consider the probability measure on Ti: 

^ , e^'^^^dx 
Pe{dx) 



/^^ e^'^^^) dx ' 

Next Lemma 14.31 is classical; we give a proof for completeness. 

Lemma 4.3 Under the assumptions of Proposition [TTT] , the following weak convergence holds: 

Pe ^ ^0 CLS 9 ^ +00 and Fg ^ 6bi as 9 ^ — oo. 

Proof of Lemma We only prove the first limit. Indeed, the second limit can be showed similarly. We 
need to show: 

Pd{A) 6o{A) as 9 ^ +oo, for any Borel set A C Ti such that ^ dA. 

If ^ A C Ti then, by assumption c{x) < c(0) for any x G A. So A CI Lt, for some t > 0, 
where /t = {x G Ti : c(x) < c(0) — t}. By assumption c is continuous at 0, so there exists an open 
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neighborhood of 0, say Vt, such that, for all x £ Vt, c{x) > c(0) — t/2. Note that, for any 9 > 0, 

< I . — dx 



- f pec{Q)-et/2 
< ^(Ft nTi)-ie-^*/2_ 

Thus, for all t > 0, \\va.Q^j^^ Fg(It) = 0. This guarantees the claim in the case when the Borel set 
j4 C Ti does not contain 0. Suppose now £ A, then ^ Ti\A and we get Pe(^) = ^ — Pe0^i\^) 
1 as goes to infinity. □ 
We can now continue the proof of the proposition. Let c{Bi) < y < c(0). By Lemma 2.3.9(6) 
in [4], we need to show that there exists a unique solution 9y of A' (9) = y. To this end, note that 
A'{9) = jj^ c{x) Pe(dx). By assumption c is continuous at and Bi, so by Lemma [4.3l and Theorem 
5.2 in pTj it follows 

lim A'{9) = c{Bi) <y <c{0) = lim A'{9). 

oo S^+oo 

Since A' is continuous and strictly increasing, the mean value theorem implies the existence and 
uniqueness of 9y. Consider now y > c(0). Note that, for 9 > 0, A{9) < 9c{0). Therefore 

9y-A{9)>9{y-cm. 

It follows that A*{y) = +00. Similarly, for y < c{Bi), we use that, for 9 < 0, A{9) < 9c{Bi) and 
deduce A*[y) = +00. Finally we prove {Hi). We first show that 

A(0/3) < A(6'), for ah > 0. (35) 

Showing (j35p amounts to show that, for all > 0, 

/ e^^(^) dx + 2/3-3 / e^^'''^/'^ dx > ^. (36) 

By Jensen's inequality it follows that 
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(the strict inequality derives from the strict convexity of the cubic power on [0,oo), and the fact 
that c is not constant on Ti). Hence the left hand side of (|36p is larger than 

9( / e^^(^)/='dx)'-3 / e^<^)/^dx + l = ^( [ e'<^y'dx-lY( [ e^^^^^'dx + l), 

and the inequality ([36]) follows. Now, let 7 < y < c(0)/3. By Theorem II. H limn^oo Pn/n = 
lim„^oo Pn/''^ = 7 < y- Thus, by Lemma 2.2.5 in [4j we have 

A*{3y) = sup(% - A{9/3)) and T{y) = sup(% - A(^)) = rjyy - A(r?^), 

where r]y is the unique positive solution of (fTO|) . Finally, ([35]) yields: 

A*(y) = yTjy - A{7]y) < yijy - A(r/y/3) < supe>o(^y " A(0/3)) = A*(3y). 
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4.2.2 Value of the Fenchel-Legendre transforms at the extrema 

In this paragraph, for the sake of completeness, we deal with the value of A* and A* at c{Bi) and 
c(0). If c is differentiable as a function from T C C to M, we denote by grad^(c) its gradient at 
X. The following proposition holds: 

Proposition 4.4 Suppose that the assumptions of Proposition \1.4\ hold and that c is differentiable 
at and Bi. If moreover, for all uj G [— 7r/2, 7r/6], grado(c) • e*'^ < and, for all uj G [27r/3, vr], 
grad^^(c) • e*"^ > 0, then 

A*(c(i3i)) =r(c(i?i)) = A*(c(0)) =r(c(0)) = +00. 

Proof. We show the proposition only for A*(c(0)). The other three cases can be proved similarly. Us- 
ing polar coordinates, we have: 

[ e^^(-')dx = [ e^^(-'")rdrda; 

for some segment I^^ = [0,a^]. The Laplace's method (see e.g. Murray [6]) gives, for all to G 
[-vr/2,V6], 

/ e'^'-'^'^^rdr ^ as 9 ^ +oo 

Ji^ ^2|grado(c) • e*"^! 

where we write / ~ 51 if / and g are two functions such that, as x ^ +00, the ratio f{x)/g{x) 
converges to 1. We deduce that, as — > +00, 

[ e^^(-)dx~e^^(o)0-2 -—l——du;. 
Jti y-7r/2 |grado(c) • e*-^! 

Since the integral in the right hand side is a finite positive constant, we have A{6) = 9c{0) — 2 log 9 + 
o(log^), and therefore 

A*(c(0)) = snp{9c{0) - A{9)) = sup(21og^ + o(log0)) = +00. 

em em 

□ 

In the next two subsections, we solve some variational problems. We refer the reader to the 
book by Buttazzo, Giaquinta and Hildebrandt [3] for a survey on calculus of variations. 

4.3 Proof of Theorem \TM.i) 

We divide the proof of Theorem ll.3( i) in 5 steps. 

Step 1: Case y ^ (c(Si)/3, c(0)/3). We have to prove that J{y) = 00. Denote by (T) C 
Mi(T) the set of probability measures on T which are absolutely continuous with respect to £. For 
a E M°'^(T), define the measures in M{,(T): 

ai{ax) = — -- — --^ ^- — --^ ^- — --^aldx), /G|l,2,3| 

Ci{x)C2{x) + Ci{x)C3{x) + C2{X)C3{X) 

where a = (1 2 3) is the cyclic permutation. Clearly ai + 02 + = a and 

$(a) < 0(ai, as, as) < c(0)/3 (37) 
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where the strict inequahty follows by assumption ^ and the fact that a is a probability measure 
on T such that a <C The above argument shows that {a £ Mf'^(T) : <I>(a) = y} = 0, for all 
y > c(0)/3. Therefore, by Theorem l4.ir z). we have J(y) = +00 if y > c(0)/3. Using assumptions ([T]) 
and dS]), one can easily realize that, for any measure P G Mft(T), > c(i?i)/3(T) and the equality 
holds only if P = 6bi- By Lemma [T2r i) we deduce that, for all a G Mi(T), 3$ (a) > c{Bi). This 
gives J{y) = 00 for all y < c(i?i)/3, and concludes the proof of this step. 

Step 2: the set function 1^ and an alternative expression for A*(3y). For the remainder 
of the proof we fix y G (c(i?i)/3, c(0)/3). For this we shall often omit the dependence on y of the 
quantities under consideration. In this step we give an alternative expression for A*{3y) that will 
be used later on. Let S C T be a Borel set with positive Lebesgue measure. Define the function of 
(?70,r/i) gR2: 

m{B,r]o,r]i) = [ e-^-'^''-'^'"^''^ dx. 
Jb 

It turns out that m{B,-) is strictly convex on (the second derivatives with respect to r/o and 771 
are strictly bigger than zero). Define the strictly concave function 

F{B, r/o, ??i) = -Vo - 3yVi - ^m{B, ??i) 

and the set function 

v{B)= sup F{B,rio,r]i). 

(r?o,r?i)eIR2 

Arguing as in the proof of Lemma 2.2.31(6) in we have: 

grad(^Q ,^^)(3m(5,-)) = (-1, -3y) u{B) = (70,71) • (-l,-3y) - 3m(S,7o,7i) 

where • denotes the scalar product on M^. Therefore, if there exist 70 = 70 (-B) and 71 = 71 (-B) such 
that 

I g-7ic(x)^^ ^ gi+70/3 and / c(x)e-^i^(^)da; = ye^+To (38) 
Jb Jb 

then it is easily seen that 

z/(S) = -(l+7o(S))-3y7i(B). 
In particular, by Proposition II. 4l fi). setting 71 (Ti) = —O^y and 7o(Ti) = h-iO^y) — 1, one has 

A*(3y) = v{Ti) = -(1 + 7o(Ti)) - 3y7i(Ti), (39) 

and 7o(Ti) and 71 (Ti) are the unique solutions of the equations in (|38p with B = Ti. Note also 
that, for Borel sets A and B such that A C C T, we have for all r/o, ?7i G M, 

m{B, m,m) - ni{A, r,o,m) = I {"^Bix) - lA(x))e-i-''«-''i^(") dx > 0. 

Jt 

In particular, for all r/o,?7i G M, F{A,r]Q,r]i) > F{B,rio,r]i). This proves that the set function is 
non- increasing (for the set inclusion). An easy consequence is the following lemma. For B CT and 
2; G C, define zB = {zx : x G B} and 

T = {Borel sets B cT: i{B) > Oand^(S n (jB)) = i{B n (fB)) = £{{jB) n {fB)) = 0}. 
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Lemma 4.5 Under the foregoing assumptions and notation, it holds: 



M{u{B) : 5 e 1} = mi{iy{B) : B £ 7 and 1{B) = 1/3} < +oo. 

Proof of Lemma The monotonicity of v implies J^(T) < i^(Ti). So the finiteness of the infimum 
follows by v{Ti) < +oo that we proved above. Note that if B e 7 then B U (jB) U (fB) C T and 
1>£{BU (jB) U ifB)) = i{B) + e{jB) + e{fB) = U{B). So 

mi{v{B) : B £7} = inf{i/(5) : 5 G T and i{B) < 1/3}. 

Now, if S G 0" is such that i{B) < 1/3, define the set C = T \ U (jB) U (fB)), note that 
£(C) = 1 - 3e{B) > and C = jC = fC. Set Ci = Cn Ti and define D = BuCi. Clearly, B C D 
and therefore > i^{D). Moreover, it is easily checked that D £ 7. Indeed, £{D) > i{B) > 

and, for instance, 

e{D n (jD)) = 1{{B u Ci) n {{jB) u (jCi)) 

< i{B n (jB)) + £(s n (iCi)) + ^(Ci n {jB)) + ^(Ci n (jCi)) = o. 

The claim follows since 



1{D) = 1{B) + i{Ci) = e{B) + i{C)/3 = 1/3. 

□ 



Step 3: the related variational problem. As above, we fix ?/ E (c(i?i)/3, c(0)/3). Recall that 
H{a I £) = +00 if a is not absolutely continuous with respect to i. So, by Theorem I4.1l fi). 

J{y)= inf H{a\£). 

Define the following functional spaces: 

"B = {measurable functions defined on T with values in [0, cxd)} 

and 

23| = |(/i,/2,/3) ^ ■■ ^ {j2fi^ = ^ and cl){£f„£f„ef,) = ^e^^ + + £f^)^ 
(recall that ^/ is the measure with density /). By Lemma l2.2( i) it follows 

J{y)= inf h(y,Mx)] (40) 

where 

= {(/i,/2,/3) e S| : m„£f„£f,)=y} 

(note that the upper script "3" in and is to remind that these spaces are defined on triplets 
of functions in 23; it is not related to the Cartesian product of three spaces). Computing the value 
of J{y) from ()40p is far from obvious, indeed Jl^ is not a convex set and the standard machinery of 
calculus of variations cannot be applied directly. The key idea is the following: consider the same 
minimization problem on a larger convex space, defined by linear constraints; compute the solution 
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of this simplified variational problem; show that this solution is in To this end, note that, again 
by LemmalZli), if {hj2,h) G % then if,{ci) = ^/.(cs) = IfM)- Therefore, we have 3i| C 8^ 
where 



= /2, /a) G : £ ( ^ /; I = 1 and, for all I G {1, 2,3},^/, (q) = y 

It follows that 



J{y)> inf H[Y,hiA 



Step 4: the simplified variational problem. Recall that y £ (c(i?i)/3, c(0)/3) is fixed in this 
part of the proof. In this step, we prove that 



I{y):= inf ^.^H\yji{x) ) (41) 
(/i>/ 

is equal to A*(3y). Clearly, the set is convex. Therefore, if is not empty, due to the strict 
convexity of the relative entropy, the solution of the variational problem (14ip . say f* = (/j^, /I , /a ) G 
S^, is unique, up to functions which are null ^-almost everywhere (a. e. ). We now compute f* and 
check retrospectively that is not empty. Consider the Lagrangian L defined by 



^(/i, /2, /a, Ao, Ai, A2, A3)(x) = /Kx) j log //(x) j + Aq h{x) - 1 

3 

+ Y,Hci{x)fi{x)-y) 



where the Aj's (i = 0, • • • ,3) are the Lagrange multipliers. For I G {1, 2, 3}, define the Borel sets: 

Ai = {xeT: fl{x) > 0}. 
Since f* is the solution of (f4T]) . by the Euler equations we have, for / G {1, 2, 3}, 

(1)1........ 

We deduce that, for all x £ Ai, 

n{x) + + n{x) = e-i-^o-A,c,(x)_ (42) 

Define the functions gi{x) := /|(ja;), (?2(x) := f^{jx) and g^{x) := fi{jx). By a change of variable, 
it is straightforward to check that {91,02,93) G and 

The uniqueness of the solution implies that a. e. 

nUx) = fiix), /Kjx) = /2*(x) and flijx) = f^{x). 
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In particular, up to a null measure set, Ai = j^^^Ai. Moreover, on ^i, the equality, a. e. Yli=i 9i{^) — 
Y^i=i fii^) applied to (|12]) gives, a. e. on Ai, exp(— 1 — Aq — A2C2(jx)) = exp(— 1 — Aq — AiCi(x)) 
(indeed x G Ai implies jx € ^2)- We deduce that A2 = Ai. The same argument on A^ carries over 
by symmetry, so finally Ai = A2 = A3. We now use the following lemma that will be proved at the 
end of the step. 

Lemma 4.6 Under the foregoing assumptions and notation, up to a Borel set of null Lebesgue 
measure it holds A\ C Ti. 

By Lemma l46l and the a. e. equality Ai = f~^A\, we deduce that A\ G T, up to a Borel set of null 
Lebesgue measure. So, by equation ([12]) and the equality Ai = A2 = A3, it follows 

fl{x) = e-^-^"-^i^(^)l(a; G Ai), a.e. 

and /|(a;) = f^{x) = f^{jx). Note that the constraints 



1 and 



read respectively 



f e-^-^«~^i^(^)d2; = 1/3 and [ c(x)e-i"^o"^i^(^) drr = y. 
JAi Jai 

This implies that the Lagrange multipliers Aq and Ai are solutions of the equations in (j38p with 
B = Ai. Moreover 

X(E^'*(^))l°g(E^'*(^)) = 3^ (-l-Ao-Aic(x))e-i-^°-^^'=(-)dx 

= -(l + Ao)-3yAi. 

Therefore (see the beginning of step 2) 



/(|:/,-,.)).o.(|:/;(.)) 



Hy) = 1 1 yjiix) I log ( > jn^) 1 dx = u{A,). 

Since yli S T we deduce that 

I{y) > mf{u{B) : B€7}. 

For the reverse inequality, take B £ 7 such that z^(-B) = sup (^^j ,j^)gR2 -F(i?, r/o, is finite. Since 
the function {r]Q,r]i) 1— > F{B,r]Q,rji) is finite and strictly concave, it admits a unique point of 
maximum. Arguing exactly as at the beginning of step 2, we have that the point of maximum is 
(70(5), 71 (-B)), whose components are solutions of equations in ([55]) . and 

u{B) = -{l + joiB))-3yji{B). 

For / G {1, 2, 3}, define the functions on T: 

gi,B ■■ X ^ e-^-^o(^)-^i(-^)'='(^)l(x G j^-^B). 
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Since 7o(-B) and ji{B) solve the equations in ([55]) . it follows easily that {gi^B^ 92,3, 93,b) £ S^- There- 
fore 



> inf Hiy^fiix)]. 



Thus 

I{y) = u{A,) = inf{z.(S) : Be7}. 

Since j4i G T, by Lemma H3] we get that (i{Ai) = 1/3. So, by Lemma we deduce that Ai = Ti 
up to a Borel set of null Lebesgue measure. Then by equation (f39l) we conclude 



I{y) = A*{3y). 

Proof of Lemma \4.6[ The argument is by contradiction. Define the Borel set 

C := {Ai n T^) U {jAi n T^) U ifAi D T^) 

and assume that £{Ai n Tf) >^0. For / e {1,2,3}, define Ai = {Ai\C) U (C n T;) and gi{x) = 
{fi{x) + + /I (x))l(x G Ai). Since Ai = f~^Ai up to a Borel set of null Lebesgue measure, 

then j^~^C = C and Ai = j^~^Ai up to a Borel set of null Lebesgue measure. So by (P2]) it follows 
that ig^{ci) = ig2{c2) = ^93(^3), and therefore 



3 / ci{x)gi{x) dx = ^ ^ (^5^ /r(^)^ dx. (43) 



Now, note that Ai C and, up to a Borel set of null Lebesgue measure, 

U ^2 U ^3 = U ^2 U yls. (44) 

So by assumption ([1]), a. e. 

3 



l(x E Ai)ci{x) < ^ l(x G A„)c™(x) 



m=l 



and the inequality is strict if x is in Cn Tz. Indeed if x G Cn Tz, then a.e. x G Am for some m ^l, 

and so q(x) < Cm(a:) by ([H). Therefore, since i.{Ai n T^^) > then i.{Cf^ Ti) > and, using (|i3|) . 
we get 

Q(a:)5/(a;) dx < ci{x)fi{x) dx = y. 



For p £ [0, 1], define the functions 

where cr = (1 2 3) is the cyclic permutation. By assumption ([6]) it follows 

ciix)gi^i{x) dx > c(0)/3 > y. 
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We have already checked that ig^ ^(q) < y, thus, by the mean value theorem, there exists p G (0, 1) 
such that (gi,p,g2,p,g3,p) G S^- The convexity of the relative entropy gives 

H{gi,p + g2,p + gi,p I < pH(gi + 52 + 53 K) + (1 - P)H{£ \ I) = pH{fl + /| + /| | (.), 

where the latter equality follows by (j44p and the definition of gi . This leads to a contradiction since 
f = {ff, /I , /s ) minimizes the relative entropy on S|. □ 



Step 5: end of the proof. It remains to check that f* = (/f , /|, /|) E Jl^. For this we need to 
prove that = f* , i f* J f*) = V- Since f* G 8^ then £j*(ci) = £/|(c2) = %(c3) = y; 

moreover, by the properties of the functions /;* it holds ^/*(q) = J^^ ci{x)fi{x) dx. So the claim 
follows if we check that 

^ef*+f*+f*)> [ ci(x)/i(x)dx. 

By Lemma [2^21 ^) we have that there exists ((71,52,53) G such that: if*^f*^f* = +^92 +^33) 
^(^/r+/2*+/3*) = <^(^<?i'^<;2>^93) and £gi(ci) = ^32(^2) = ^93(c3)- In particular, 

3 3 /i 3 

3^(^/r+/2*+/3*) = X] / c/(2;)5Ka:^)dx = ^ / ^ci{x)gi{x)dx 

3 „ 3 

Crn(2;)^5/(2;)dx (45) 

m=l •''^^ 1=1 
3 

<^m{x)f^{x)dx 

m=l -^^^ 

ci(x)/i (x) dx 

where in ()45p we used assumption ([T|). This concludes the proof of Theorem II. 3l fi). 



4.4 Proof of Theorem \1.3{ n) 

Some ideas in the following proof of Theorem ll.3( ii) are similar to those one in the proof of Theorem 
ll.3l fi). Therefore, we shall omit some details. We divide the proof of Theorem ll.sr ii) in 3 steps. 



Step 1: Case y ^ (c(i?i)/3, c(0)). As noticed in step 1 of the proof of Theorem ILSf i). for any 
measure [3 G Mfe(T), /3(q) > c(i?i)/3(T) and the equality holds only if /3 = 5bi - We deduce that, for 
ah a G Mi(T), 3^'(a) > c{Bi). Therefore, by Theorem [iT^ii) , J{y) = +00 if y < c(Bi)/3. Now, 
note that, for a G Mi(T) it holds 

^{a) = max q(x)a(dx)^ < c(0) max a(T;) < c(0) 

where the strict inequality follows by assumption ^ and a <^ i. Therefore, using again Theorem 
l4.1l fM). we easily deduce that J{y) = +00 if y > c(0). 
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Step 2: the set function ^u. For the remainder of the proof we fix y S (c(i?i)/3, c(0)), and 
we shall often omit the dependence on y of the quantities under consideration. In the following we 
argue as in step 2 of the proof of Theorem ll.3l fi). Let i? C T be a Borel set with positive Lebesgue 
measure and define the function of (r/o -.Vi) 

q{B, r/o, m) = 2e-^-''^H{B n T2) + f Q-i^m-mc{x) 

JBnTi 

Clearly, q(B,-) is strictly convex on M?. Define the strictly concave function 

G{B, r]o,r]i) = - ym - q{B, r/i) 

and the set function 

H{B) = sup G{B,r]o,r]i). 

(»?o,m)eK2 

If there exist 79 = 7o(^) and 7^ = 7i(-B) such that 

[ e-^^^^^^dx + 2£{B D T2) = e^+^o and / c(x)e-^i"(^)dx = ye^+~o (45) 
JBnTi JBnTi 

then we have 

^{B) = -{l+%{B))-y^,{B). 
In particular, by Proposition II. 4( u). setting 71 (T) = —r]y and 70 (T) = A(r/y) — 1 one has 

r(y)=/i(T) = -(l + 7o(T))-y7i(T) if 7 < y < c(0) (47) 

and 7o(Tr) and 7i(T) are the unique solutions of the equations in (j46p with B = T. Recall also that 
in step 2 of the proof of Theorem II. 3( i) we showed: 

A*{3y) = -{!+ 7o(Tri)) - 3y7i(Ti) if c(Bi)/3 < y < 7 

where 70 (Ti) and 71 (T) are the unique solutions of the equations in ([38]) with B = Ti. Note 
that, for Borel sets A and B such that yl C C T, we have, for all r]Q,r]i £ M, G{A,r]Q,r]i) > 
G{B,rjQ,r]i). This proves that the set function fi is non-increasing (for the set inclusion). An easy 
consequence is the following lemma: 

Lemma 4.7 Under the foregoing assumptions and notation, it holds: 

mi{fi{B) : BCT}=T{y) if ^ < y < c(0) 

Step 3: the related variational problem. As above we fix y G (c(-Bi)/3, c(0)); as in the proof 
of Theorem II. 3r i) we denote by H the set of Borel functions defined on T with values in [0, 00). By 
Theorem I4.1( ii). we have 



m = -f^^H{f) 



where 



U = |/ E B : ^(/) = 1 and ^^f^^i^j ci{x) f {x)dj^ = y^ . 
Note that / € IX if and only if the functions x i— > f{jx) and x ^ fij'^x) are also in U and so 



J{y) = MH{f) (48) 
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where 

V = {/ G 03 : l{f) = 1, if^j^ici) = y, £f^j^{c2) < y, if^j^ic^) < y} 



The optimization problem (j48p is a minimization of a convex function on a convex set defined by 
hnear constraints. Thus it can be solved explicitly. Therefore, if V is not empty, since the relative 
entropy is strictly convex, the solution of the variational problem ([18|). say /* G V, is unique, up to 
functions which are null ^-almost everywhere. We will compute /* and show that V is not empty 
at the same time. So assume that V is not empty and define the function 

g{x) = r(x)iTi(x) + r{jx)iT,{x) + nfx)ij,ix). 

It is easily checked that g gV and H{g) = H{f). The uniqueness of /* implies that 

for almost ah x G T2, = f*{x). (49) 

Therefore, up to modifying /* on a set of null measure, /* G V where 

V = {/ G 03 : l{f) = 1, ^/nr^(ci) = y, if^j^{c2) < y} 

and the variational problem reduces to J{y) = infjgy/ H{f). Consider the Lagrangian L defined by 

Hif, Ao, Ai, A2)(x) =/(x) log fix) + Ao(/(x) - 1) 

+ Ai(ci(x)/(x)1ti(x) -y) + X2{c2{x)f{x)lj^{x) - y) 

with 

I C2{x)r{x)dx-^ =0. 



The two cases A2 = (i. e. /* is not constrained on T2) and A2 7^ (i. e. /* is constrained on 
T2) are treated separately. For each case, we solve the variational problem. The optimal function 
is denoted by fu for A2 = and by /c for A2 / 0, so that /* = avg mm{H{ f u), H{fc)). Assume first 
that A2 = so that /* = fu and define the Borel set: 

Au = {x£T: fu{x) > 0}. 

By the Euler equations we get, for all x G T, 

= lT,nA„(x)e"i-^«-^i^^(=^) + l(Tr,uT3)nA.(^)e-^-^°. (50) 

By (fl9l) we have i{Au n T2) = i{Au n T3), and so the constraints i{fu) = 1 and ^/uj-j^ (ci) = y read, 
respectively, 

/ e~^i"(^) dx + 2i{Au n T2) = e^+^« and / c(x)e-^i^(^) dx = ye^+^° . 

With the notation of step 2, this implies that Aq = Jq{Au) and Ai = ji{Au) are the solution of the 
equations in (j46p with B = A^. In particular, 

^i{Au) = -(1 + 7o(^«)) - yiMu) = H{fu) 

where the latter equality follows from the computation of the entropy using ()50p . By Lemma 14.71 
we deduce that 

H{fu)>T{y) if7<y<c(0). 
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By (|Tr|) we have H{h) = A (y), where 

h{x) = lj,{x)e-^-^'>m-^imc(x) ^ ^^^^^^^^^g-i-7o(T) 

and 7o(T), 7i(T) are the unique solutions of the equations in (|46p with B = T. Now we prove that 
/i G V, for 7 < y < c(0), so that 

H{fu)=Tiy) if7<y<c(0). (51) 

Recall that — 7x(T) is the unique solution of 

/tt^ c(a;)e^^(^) dx 
/^^ e^^(^) dx + 2/3 ~ ^' 

The function 

L c(x)e^^^^^ dx 
e^'=(^') dx + 2/3 

is strictly increasing (as can be checked by a straightforward computation) and, for ^ = 0, it is 
equal to 7. Therefore, since 7/ > 7, we have — 7x(T) > 0. It implies that 

/ c(x)e-i-7om-7i(TMx)^^^y > f c(x)e-i"^omdx = 7e-i-7om. 
In particular, /i G V. Now we deal with the case A2 / 0. We have 

In particular, if we set fc,i{x) = l(x G T/)/c(x), we get (/c,i, /c,2, /c.a) G §0- By step 4 of the proof 
of Theorem II .3^ 7) ■ it implies that 

H{f,) > inf H{h + h + h) = A*(3y) = Hif^ + Z^* + f^), 

(/l,/2,/3)eS3 

where f* = (/ij/^j/s) was defined above. Since /i + /I + /a G '^^j we deduce directly that a.e. 
fc = ft + /I + /s* and 

H{f,)=A*{3y). (52) 

It remains to find out for which values of y the Lagrange multiplier A2 is equal to zero. First of all 
note that if y = 7 then the function identically equal to 1 is in V. We deduce that /* = 1 and so 
A2 = (since the optimal solution is not constrained on T2) and J{'y) = = A* (37). Now assume 
7 < y < c(0). By Proposition [Olf ziil . we deduce A*(y) < A*(3y). It follows by ^B) and ([52D that 
H{fu) < H{fc). Recall that /* = arg min(i/(/u), i7(/c)), thus A2 = and J(y) = A*(y). It remains 
to deal with the case c(i?i)/3 < y < 7. The following lemma holds: 

Lemma 4.8 Under the foregoing assumptions and notation, if c{Bi)/3 < y < 7 then J{y) > J{y)- 
Then, by Theorem II. 3( i) and ()52p we get 

A*(3y) = J(y) < J(y) = mm{H{f^), H{f,)) < A*(3y). 
This concludes the proof. 
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Proof of Lemma \4 ■ 4 Choose y < z < 'j. By construction P{Pn ^ nz) < P{pn < nz). Taking the 
logarithm, applying Theorem 14.11 and recalling that J(y) = J{y) = +00 for y < c{Bi)/2> we have 

— inf Jit) < liminf — log P(p„ < nz) < limsup — logPfpr, < nz) < — inf J{t). 



Therefore 



Jiy) > inf J(t) > inf J(t) = J(z) 

iG(c(Bi)/3,^) te(c(Bi)/3,z] 



where the latter equality follows since J{y) = A*(3y) is decreasing on (c(i3i)/3, 7). Recalling that 
J(y) = A*{3y) is also continuous on (c(Si)/3,7), the claim follows letting z tend to y. □ 



5 Model Extension 

5.1 The analog 1-dimensional model 

The analog 1-dimensional model is obtained as follows. There are n objects on (0, 1), say {1, . . . , n}, 
and two bins located at and 1, respectively. The location of the k-th object is given by a r.v. Xk 
and it is assumed that the r.v.'s {Xk}i<k<n are i.i.d. and uniformly distributed on [0, 1]. The cost 
to allocate an object at x S [0, 1] to the bin at 0, respectively at 1, is c{x), respectively c{l — x). The 
asymptotic analysis of allocations which realize the optimal and the suboptimal load can be carried 
on using the ideas and the techniques developed in this paper. Due to the simpler geometry of 
the 1-dimensional model, many technical difficulties met in the 2-dimensional case disappear, and 
with the proper assumptions on the cost function, it is possible to state and prove the analog of 
Theorems [Til O and \TM 



5.2 Random cost function 

An interesting and natural extension of the model takes into account random cost functions. Let 
2, be a Polish space and = {k = l,...,n) a r.v. taking values on Z^. Assume 

that: the sequences {^fc}i<fe<n and {Zfc}i<fc<„ are independent; the r.v.'s {Zfc}i</;<„ are i.i.d. 
with common distribution Q; the r.v.'s Z^, Zf and Zf are i.i.d.. Let c : T x 2.^ — > [0, 00) be a 
measurable function. We consider an extension of the basic model where the cost to allocate the 
fe-th object to the bin at Bi [l = 1,2,3) is equal to q(Xfc,Zfc). Here, for z = {z^ , z^ , z^) , the cost 
functions are defined in such a way that they preserve the spatial symmetry: ci(a;,z) = c(x,z), 
C2(x,z) = c(j^2;, (z^, z^, z^)) and C3(x,z) = c{jx,{z'^ , z^ , z"^)). The load associated to an allocation 
matrix A G An is 

Pn{A) = max i^^Y^auci{Xk,7jk)^ . 

In a wireless communication scenario we have 2. = M+, and the typical cost function is of the form 

a -|- min{6, — i?2|~"} + niin{6, z^jx — -Bsl^"^} 

where a > 0, a > 2 and b > (A\/3/2)~°. The additional randomness in the cost function models the 
fading along the channel (see e.g. [9]). The suboptimal allocation A = {ak^i)i<k<n,i<i<3 is obtained 
by allocating each point to its less costly bin. To be more precise, assume that i (8) Q-a.s., for any 
/ / m, q(x,z) 7^ Cm{x,z). Then, setting 

ak,i = l{ci{Xk,Zk) < m.inCm{Xk,Zk)), 
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the suboptimal allocation matrix is a.s. well-defined. Consider the suboptimal load p„ = Pni^) 
and the optimal load pn = min^gyi^ Pn(^)- Exactly as in the proof of Theorem 1 1.1^ one can prove 
that, a.s. 

,in, ^ii = 5. = / ,(c,(..z)<,nin.„.(.,z))d.«dz). 

n^oo n n—too n JjxZ^ ™7^' 

Deriving analogs of Theorem 11.21 and Theorem 11.31 is an interesting issue. For the central limit 
theorem, an analog of the suboptimal allocation matrix A in Proposition 13.11 should be defined. For 
the large deviation principles, the contraction principle can be applied as well, but it might be more 
difficult to solve the associated variational problems. 

5.3 Asymmetric models 

Most techniques of the present paper collapse when the symmetry of the model fails, e.g. the region 
is not an equilateral triangle, the locations are not uniformly distributed on the triangle, the cost 
of an allocation is not properly balanced among the bins. For a result on the law of large numbers 
in the case of an asymmetric model, we refer the reader to Bordenave [2]. 

6 Appendix 

6.1 Proof of Lemma 12.11 
Continuity of (p. By the inequality 

I max{ai,a2, 03} -max{6i, 62, h}\ < |ai - ^i| + 1^2 - &2I + - ^sl, for all 01,02,03, 61, 62, ^3 > 0, 
we get 

10(01,02, as) - 0(/3i,/52,/?3)| < |ai(ci) - /3i(ci)| + |a2(c2) - (32{c2)\ + {asics) - /33(c3)|. (53) 

Since c is continuous, if the sequence {{ai,a2,a^))n>i G 3V[;,(T)^ converges to (/3i, /32; /^s) (with 
respect to the product weak topology), then 

lim |a^(ci) - Pi{ci)\ = lim |a^(c2) - P2{c2)\ = 

n— >oo n— KX) 

and 

lim |a^(c3)-/33(c3)| =0. 

n^oo 

The conclusion follows combining these latter three limits with (j53p . 

Continuity of ^. For each / E {1,2,3}, the projection mapping a 1— > a^j^ is continuous. Hence, 
the continuity of ^ follows by the continuity of (p. 

Continuity of $. Note that, for each fixed a G Mi(T), it holds 

<I>(q) = (/)(ai, 02; as) for some ai, 02, as G Mfe(T) : ai + 02 + = a 

(indeed, the set {(ai,a2,a3) G Mb(T)^ : ai + a2 + as = a} is compact with respect to the 
product weak topology and the functional (p is continuous). For each integer K > 0, consider the 
open covering of T given by the family formed by the open balls centered at x G T with radius 
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XjK. Then by a classical result (see, for instance, Proposition 16 p. 200 in Royden [8j)) there exists 
a finite collection {V'ra}i<n<7V of continuous functions from T to T such that 



N 



tpnix) = 1 for each x G T, £{supp{ipn)) < 2/ir for each n = 1, . . . , N. 



n=l 



Here the symbol supp{ipn) denotes the support of Vn- Let / be a continuous function on T, 
consider the modulus of continuity of / defined by ws{f) = sup|5_t|<5 — f{t)\, and set 
/„ = sup^.g^yppj-^^) f{x). Note that, for all measures fi G 'Mb{T), 



N 



N 



^ — ' K ^ — ' K 



(54) 



For i = 1, 2, 3, define r 



n=l n=l 



if a{^l^n) > and = otherwise. Moreover, for /3 G M;,(T), set 



TV 



A(dx) = r>n(x)/3(dx), i = 1, 2, 3. 



(55) 



n=l 



Since al{^jJn) + «2(V'n) + a3(V'n) = a(V'n)! by the properties of the sequence {'4'n}i<n<N we have 
Pi + P2 + = P- For any continuous function / on T we have, for i = 1, 2, 3, 



N 



n=l 



- a^{f)\ =\ E i^W^n) - a,(/Vn)) 
n=l 

<| E< (^(^^«) - "(^^")) + I E< (M^n) - aifiJn)) 

n=l 

N 



N 



(56) 



n=l 



Note that < 1, and therefore 



AT 



Vr;(/3(/V^„) -"(/V'n)) < iV max |/3(/V^„) -a(/Vn)|. 

— ' l<n<N 



n=l 



(57) 



Using again that < 1 and ([Ml) with // = a, we have 



(58) 



n=l 

By the definition of rj^ and (j54|) it follows 

N 



n=l 



N 



J2 i'^nfna{lpn) - Ctiiflpn)) \ = \^ UnaMn) - ai{f^n)) < W2_{f). 



n=l 



n=l 



Cohecting 1^, ([57]), §^ and ([59]) we have 



mf) - «.(/)! < Nms.x IPifiJn) - a(/V'n)| + 2w^{f). 

l<n<N K 



(59) 



(60) 
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Now, let {Z?™} C Mi(T) be a sequence of probability measures converging to a for the topology of 
the weak convergence. We shall prove 

lim = $(a). 

m— >oo 

We first prove 

limsup ^>(a). (61) 

m^oo 

Let K be as above and define the Borel measure as in (j55|) . with fi^ in place of /3 (the definition 
of r\ remains unchanged). By inequality (|60p and the weak convergence of /3"^ to q, it follows 

limsup|/3r(/) -««(/)! <2«;^(/). 

m— >oo 

Applying the above inequality for / = ci , / = C2 , / = C3 and using the inequality ([53]) , we get 

limsup|(/.(/31",/32",/33") -'/'(«i>«2,a3)| < 6t(;^(c). 

m— ►oo * 

Note that by the definition of ^ and the choice of the a^'s, <I>(a) = (^(ai, «3) and <I>(/3™) < 
0(/3J",/3^,/3£^), therefore 

lim sup ^> (/?"") < $(a) +6z/;a(c). 



The above inequality holds for all and letting K tend to infinity, we obtain (16ip . We finally 
check the lower semi-continuity bound 

liminf ^>(/3'") > ^>(a). (62) 

m^oo 

Arguing as at the beginning of the proof, we have, for each fixed m > 1, 

$(/3'") = 0(/3f , /3™) for some , , /^g"^ G Mf,(T) : /?|" + /?™ + = /?"^. 
Now, consider an extracted subsequence {jnk)k>\ such that 

liminf$(/3'^) = lim (/.(/3™\ ;3^'=). 

m^oo fc^oo 

As already pointed out, Mf,(T)'^ is compact with respect to the product weak topology. Therefore, 
up to extracting a subsequence of (mfc)fc>i, we may assume that (/3^* , Z?^*" , ) converges to 
(/?i, /32, /Ss) G Mfe(T)3. By construction, + + f3f = (3"" and /?™ converges to a, thus, we have 
A + /?2 + /^a = a- Then the definition of $ gives 

HPi,P2,P3)>Ha). 

Also the continuity of (j) implies 

lim </)(/?r,/32™N/33"''=) = </)(A,/32,/33). 

fe— >oo 

The matching lower bound (j62p follows. 
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6.2 Proof of Lemma 12.21 

Proof of (i). For each a G Mi(T), the set 

{(ai, 02, as) G Mb(T)^ : ai + 02 + "3 = "} 

is convex; moreover, the functional (/> is convex on Mft(T)'^. Therefore, by a classical result of convex 
analysis, there exists, (ai,a2,a3) £ Mb(T)^, such that <I>(a) = i;^(ai, 02, 03)- 

In order to prove that ai(ci) = a2(c2) = oslcs), we reason by contradiction. Assume, for 
example, that <I>(a) = ai(ci) > max(Q;2(c2), 03(03)). For p G (0,1), define (/3i , /32 > /93 ) = (pQ;i)(l — 
p)ai + 02, 03). We have I3i + (^2 + l^i, = ot and 

</'(/3i,/?2,/?3) = max(pai(ci), (1 -p)ai(c2) + 02(02), 03(03)). 

In particular, for p large enough, (piPi^ (^2, Ps) = poi(ci) < 0(01,02,03). This is in contradiction 
with <I>(o) = i?!)(oi, 02, 03). Now, assume, for example, that <I>(o) = oi(ci) = 02(02) > 03(03). The 
same argument carries over, by considering, for p G (0, 1), (/3i, /52, /93) = (pai,po2, 03 + (1 —p){ai + 
03)). All the remaining cases can be proved similarly. 

Proof of (ii). Since An C S„, we have pn < Pn, and therefore we only need to establish the 
claimed lower bound on Let B* be an optimal allocation matrix for pn and define the set 

/ = {A; G {1, . . . , n} : there exists / G {1, 2, 3} such that bli G (0, 1)}. 

Define the matrix A = (a^/) G An by setting a^i = 6^^, for any / G {1, 2, 3}, if A; ^ /, and a^i = 1, 
afe2 = CLk3 = if A; G /. Letting |/| denote the cardinality of /, we have 

Pn = max I bliCiiXk) + ^ bliCi{Xk) j 

\kel kfl / 

> max y^afcic(Xfc) + ^afcic(Xfc) - |/|||o||oo, max ^afc;c/(Xfc) 

Vfce/ k0 '^^^'^Hfc^/ // 

n \I\ ||c||oo- 



\k=l 



Thus, the claim follows if we prove that |/| < 3. Reasoning by contradiction, assume that |/| > 4 
and, for j = 1, 2, 3, 4, denote by kj G / four distinct indices in /. For each kj there exists Ij G {1, 2, 3} 
such that b*j^,i, G (0,1). Since 

b* + b* =1 

mg{l,2,3}\M 

we deduce that there exist rrij G {1,2,3} \ {Ij} such that b*^^^_ G (0,1). Thus if |/| > 4, there 
exist distinct ki,kj G {l,...,n}, distinct Zi,mi G {1,2,3} and distinct lj-,fnj G {1,2,3} such that 
bk,h-,bk,m,-,bkjij-,bkjni^ G (0,1). Choose e G (0, min{5^^;^, 6^ imi^b%.i-,b*^_^^}) and define the matrix 
= {bli) G 03, by 
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and = 6^^ otherwise. We define similarly B ^ hy replacing e by —e. By part (i) of the lemma, 
the optimal allocation matrix B* satisfies 



- \A:=1 fc=l / fc=l 

n n 

k=l k=l 

Therefore 

= max(^e(Q,(XfcJ - q^. (X^J), ±e(c™^ (X^J - c^^{Xk^))) > 0. 

It gives q.(Xa;.) = ci^{Xkj) and Cmi(-'^A:i) = Cmjl-'^fcJ but it a.s. cannot happen since, by assump- 
tion, £{c-^{t})) = for all t > 0. 

Proof of (iii). It is an immediate consequence of (ii). 

6.3 A particular cost function: the inverse of signal to noise plus interference 
ratio 

In this subsection, we prove that the following cost function 

, , a + min{6, |x — i?2|~°l + min{^) 1^; — i?3|~"} ^ 

where a > 2, a > and 6 > (A\/3/2)^", satisfies ([I]), ([2]), ([3]), (jH) and (p. To avoid lengthy 
computations we only checked numerically the first inequality in ([6]). The typical shape of the 
function 

^/ ^ ^ Ci{x)c2{x)c3{x) 

Ci{x)c2{x) + Ci{x)c3{x) + C2{x)c3{x) 

is plotted in Figure O which shows that L attains the supremum at x = 0. Finally, we show that, 
for fixed a > 2 and a > 0, for all b large enough, the second inequality in ^ holds. 

We first check assumption ([1]). We consider only the case 1 = 2, being the case I = 3 similar. Let 
X G T be such that |x — < |x — i?2|- Then necessarily, |x — i?2| > A\/3/2. With our choice of b, 
we deduce that 

min{6, \x - SsT"} = \x - ^al"" < min{6, \x - 



By construction 



a + min{6, |x — "} + min{6, |x — -B3I "} 
mmjo, \x — B2\ j 



and so ([T]) follows easily. 

It is immediate to check that c is a Lipschitz function, and the axial symmetry around the 
straight line determined by and Bi maps B2 into B3. Thus assumptions ^ and ^ follow. 



35 



In order to check ([S]), we note that if x G Ti, then, for I = 2,3, \x — Bi\ > |x — Thus, for 
/ = 2,3, min{6, \x — Bi\~°'} < min{6, \x — Bi\~"}, and we deduce 

a + min{6, |x — i?2|~°} + niin{6, |a; — Ssl""} 

~ min{6, |x - -Bi|-°} 

a 

< l_ 2 

~ min{5, |a; — -Bi|~"} 

< A"a + 2 = c(0), 

where the last inequahty is strict if x 7^ 0. Similarly, a + min{6, \x — -62!""} + min{5, \x — B^l^"'} is 
minimized for x = Bi and min{5, |x — -Bi|~°} is maximized for x = Bi. So, for x ^ Bi, c(x) > c{Bi). 
Now we check assumption ([3]). Define 

Ai = {xGT:\x-Bi\<b-^/'^}, 1 = 1,2,3. 

With our choice of 6, if / 7^ m, we have Ai n Am = 0- Define 

^0 = T\ {A, U A2 U A3) . 

Note that, by construction, on each set Ai, / = 0, 1, 2, 3, the sign of 6 — |x — " is constant for 
each m = 1, 2, 3. To prove ([3]), we shall check that, for all t > and I = 0,1, 2, 3, 

e{Ainc-\{t})) =0. (63) 

We shall only prove the above equality for I = 0, the other cases can be shown similarly. Note that 

c{x) = \x- 5ir(a + \x- B2\~'' + \x- ^sl""), Vx G ho- 
using polar coordinates we have 

r2n roo 

£{Aonc-\{t})) = de l{re'^ £ Ao}l{c{re'^) = t}r dr. 
Jo Jo 

We shall check that, for an arbitrarily fixed 9 £ [0, 27r), the function 



is strictly monotone, where 

le = {r : r >0, re*^ G T}. 

So, for any fixed 6 G [0, 27r), the function l{re*^ G Ao}l{c(re*^) = t} is different from for at most 
one r, and therefore the equality ()63p for I = follows. In the following we shall only prove that cq 
is strictly decreasing on Iq for 6 G [— 7r/6, 7r/6], the other cases can be treated similarly. First, note 
that since 6 G [— 7r/6, 7r/6], as r increases, \re^^ — Bi\°' decreases, while |re*^ — i?3|" increases. Thus, 

r I— > a|re*^ — Bil" and r (j^mz^^ decreasing. Note also that, for 9 G [— 7r/6,0], as r 

increases, |re*^ — -B2I" increases. Thus it suffices to prove that, for a fixed G (0,7r/6], the function 



0, A (2cos ( ^ 
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is non-increasing. Consider the orthonormal basis {61,62} with ei = e*6 and 62 = e *3. Setting 
/? = 7r/6 - 6* G [0,7r/6), yi = A/2 and y2 = A^/3/2, we have 

re*^ = r cos/?ei + r sin/?e2, Bi = j/iei + 7/262, -62 = 7/161-1/262 

and 

^ _ [yi -rcosP f + (j/2 - r sin P f 
(2/1 — r cos + (2/2 + sin /3)2 

The derivative L'g{r) of Lg{r) has the same sign of 

— (cos P{yi — r cos j3) + sin/3(y2 — r sin/?)) ((yi — r cos /3)^ + (?/2 + sin /?)^) 
+ (cos f3{yi — r cos f3) — sin/3(y2 + r sin/3)) ((yi — r cos /?)^ + {y2 — r sin/3)^) . 

After simphfication, we get easily that L'g{r) has the same sign of 

— 2r cos /? sin /3 — ((yi — r cos /3)^ + 2/2 ~ sin^ /3) sin (3. 

This last expression is less than or equal to 0. Indeed, for r G [0, A(2 cos /3)~"^], we have < rsin/3 < 
2/2- Hence Lg is non- increasing on its domain. 

Finally, we check that, for fixed a > 2 and a > 0, it is possible to determine b > (A\/3/2)~" so 
that the second inequality in ^ holds. We deduce 



/ 



, a + min{6, |x-52| °} + min{6, |x - B3I "} 

c(2;)dx> / r= dx (64) 

~ k (Ay3/2)- 

a + min{6, \x - ^2|~"} + \x - ^al"" ^g^^ 



2 



(AV3/2) 



> f +-^"/^"^/%(AV3/2)- / |x-i.3|-"dx. (66) 
- (AV3/2)- (A^/3/2)- ^ ^ Wt2 

Here ([Ml) and ([65]) follow since on T2 we have \x - Bi\-'^ < (A\/3/2)-° < 6 for Z = 1,3; ([66]) is 
consequence of the inequality \x — i?2|~° > b, for any x G A2 H T2. The claim follows noticing that, 
due to our choice of a, c(0)/3 is strictly less than the quantity in ([66j) . for b large enough. 
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Figure 3: The function L with a = 2.5, a = 1 and b = 10. 
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